Pair production in laser fields oscillating in space and time 
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The production of electron-positron pairs from vacuum by counterpropagating laser beams of 
linear polarization is calculated. In contrast to the usual approximate approach, the spatial depen- 
dence and magnetic component of the laser field are taken into account. We show that the latter 
strongly affects the creation process at high laser frequency: the production probability is reduced, 
the kinematics is fundamentally modified, the resonant Rabi-oscillation pattern is distorted and the 
resonance positions are shifted, multiplied and split. 
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In the presence of very strong electromagnetic fields the 
quantum electrodynamic vacuum becomes unstable and 
decays into electron-positron (e + e~) pairs This phe- 
nomenon has been realized experimentally, e.g., in rela- 
tivists heavy- ion collisions Q. Shortly after the inven- 
tion of the laser almost 50 years ago, theoreticians began 
to study e + e~ pair production by intense laser light [3[. 
Because of the remarkable progress in laser technology 
during recent years, the interest in the process has been 
revived since an experimental investigation of pair cre- 
ation by pure laser light is coming into reach. The only 
observation of laser-induced pair production until now 
was accomplished ten years ago at SLAC (Stanford, Cal- 
ifornia) , where a 46 Ge V electron beam was brought into 
collision with an intense optical laser pulse In this 
experiment, a 7-photon produced via Compton scatter- 
ing or the electron Coulomb field assisted the laser beam 
in the pair production. 

The most simple field configuration for realization of 
purely laser-induced pair production consists of two 
counterpropagating laser pulses of equal frequency and 
intensity (see [5|, |y, |7|, lS( and references therein). The 
resulting standing wave is inhomogeneous both in time 
and in space which represents a formidable task for the 
nonperturbative quantum field theory, see e.g. [9]. All 
theoretical investigations so far have approximated the 
standing laser wave by a spatially homogeneous electric 
field oscillating in time. This dipole approximation is ex- 
pected to be well-justified in optical laser fields, where 
the wavelength is much larger than the typical length 
scale of the process: A ^> I ~ 2m/(\e\E) in natural 
units (h = c = 1) which are employed throughout. In 
terms of the relativistic parameter £ = \e\E/(mw), this 
relation corresponds to £ ^> 1. Here E and uo are the 
laser field and its frequency, e and m the electron charge 
and mass, respectively. Meanwhile the experimental re- 
alization of laser-induced pair production is also exten- 
sively discussed in connection with upcoming x-ray free- 
electron laser (XFEL) facilities @, Hj]. In this case, how- 
ever, the laser frequency is high, £ < 1 and the magnetic 
field component is not negligible. The latter, in general, 
can have an important influence on the pair creation pro- 



cess. This is most evidently demonstrated by the fact 
that a single plane laser wave cannot extract pairs from 
the vacuum, whereas a purely electric field can. 
In this Letter, we calculate e + e~ pair creation by two 
counterpropagating strong laser pulses (CLP) of high fre- 
quency (uj < m and £ < 1), taking into account explic- 
itly the temporal and the spatial dependence of the laser 
fields and their magnetic components. The pair produc- 
tion is modelled via an electron transition from an initial 
negative-energy state to a final positive one, based on a 
solution of the Dirac equation [10[ . By numerically solv- 
ing the latter, we demonstrate that the laser magnetic 
field component significantly reduces the pair creation 
rate. Moreover, it leads to distinct qualitative changes 
in the process as the resonant Rabi oscillation pattern is 
strongly disturbed. These features can be used to test 
strong-field QED in temporally and spatially inhomoge- 
neous fields. 

We employ an advanced computer code that solves the 
Dirac equation in an arbitrary external potential on a 
two-dimensional spatial grid [jjj . An initial wave-packet 
ty(t = 0) in the negative-energy continuum, representing 
an electron in the Dirac sea, is propagated via the split- 
operator algorithm. Under the influence of the external 
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FIG. 1: Pair creation in an oscillating electric field with fre- 
quency u = ra/200 and critical field strength E c = m 2 /e. 
Three snapshots of the probability distribution are taken at 
times as indicated from the evolution of an initially negative- 
energy Gaussian wave-packet at rest at the origin. The length 
is scaled in multiples of the Compton wavelength A c = 1/ra. 
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FIG. 2: (color online) Disturbed Rabi oscillations: Pair pro- 
duction probability versus pulse length for fixed turn-on and 
turn-off phases of sin 2 -shape and a duration of half a cycle 
each at £ = 1. (a) The red dashed line and the black solid 
line show the OEF case for uo — 0.49072m and the CLP case 
for uo — 0.4721m, respectively, each corresponding to a 5- 
photon resonance, (b) CLP case for various frequencies. The 
black solid curve coincides with the resonant case in part (a), 
while the others correspond to off-resonant situations. 



FIG. 3: (color online) Resonant probability spectrum: Maxi- 
mal value of the pair production probability during Rabi os- 
cillation at £ = 1, varying the pulse length up to 200 cycles. 
The red crosses show the OEF spectrum; the peak labels de- 
note the absorbed photon number. The black triangles show 
the CLP spectrum. Here the labeling signifies the number 
of absorbed photons from the right- left propagating waves. A 
splitting occurs, as indicated by arrows for the example of the 
(3-2) peak. The frequency axis is plotted reciprocally [l2[. 



field an e + e~ pair can be produced. The transition am- 
plitude is determined by projection of the wave function 
onto positive-energy states 3>p1~^ after the external field 
has been turned off. In this approach an intuitive graph- 
ical interpretation for the production process is possible. 
As an example Fig. [1] shows the time evolution of an ini- 
tially negative-energy Gaussian wave packet at rest under 
the influence of an oscillating electric field (OEF). When 
the e + e~ pair is produced, a droplet is separated from 
the wave packet which moves opposite to the initial one. 
The droplet is a positive energy state and represents the 
created electron. The change of the sign of energy is evi- 
denced by the change of the group velocity of the droplet 
wave packet. 

Our numerical approach enables us to take into account 
the full space-time dependence of the laser potential 
and to investigate the effect of the magnetic field. Be- 
fore elaborating on the differences of the pair produc- 
tion dynamics in an OEF and CLP, we briefly summa- 
rize the main results for pair production in an OEF. 
Due to momentum conservation in OEF, the problem is 
reduced to a two-level system consisting of a negative 
and a positive-energy state where a multiphoton reso- 
nance 2<2o(p) = nuo enforced by energy conservation can 
exist, leading to Rabi oscillations between these states 

@. Here q (p) = (1/T) J Q T dtyj (p - eA(t)f + m 2 is the 
quasi-energy of the laser-dressed state, p the canonical 
momentum, T the pulse duration and n the number of 
absorbed photons. The relativist ic parameter £ distin- 
guishes between different interaction regimes of pair cre- 
ation in laser fields. For ^ « 1 the process probability 
follows a perturbative power law dependence: W ~ £ 2n , 



whereas for £ ^> 1 (E < E c ) a tunneling behavior: 
W ~ exp(— ttE c /E), where E c = m 2 /\e\ is the critical 
field Most interesting is the nonperturbative regime 
(£ < 1), where non-dipole dynamics is pronounced and 
no simple asymptotic formulas exist. Therefore, we re- 
strict the following discussion to £ = 1 (for both laser 
beams taken together). 

Now we turn to pair production in CLP. Our calculation 
always assumes an initially narrow Gaussian wave-packet 
centered around p = of width Ap « am <C m, with 
the fine structure constant a, lying in the negative-energy 
continuum. A particular spin state along the magnetic 
field direction is chosen; the opposite spin orientation 
would give identical results. The wave-packet is exposed 
to two linearly polarized laser pulses counterpropagating 
along the z-axis and featuring sin 2 turn-on and turn-off 
phases of half a cycle each. The number of plateau cy- 
cles with constant intensity is variable. To save com- 
puting time the fields are switched off when they do not 
overlap anymore. The pair production probability is ex- 
pected to depend on the pulse length as known from the 
OEF case. The latter is shown in Fig. 0(a), where the 
red dashed curve shows the expected Rabi oscillation in 
an OEF for uo = 0.49072m corresponding to an n = 5 
photon resonance, because the value of the quasi-energy 
is go(0) ~ 1.21m which is independent of frequency for 
given £. To compare this with the CLP case one has to 
find the n = 5 resonance frequency. In this case qo is 
no longer analytically computable. Fig. [2](b) shows the 
result for various frequencies starting from uo = 0.46m on 
a logarithmic scale. The resonance frequency for n = 5 
turns out to be uo = 0.4721m. The identification of res- 
onances in the CLP case will be explained below. The 
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FIG. 4: (color online) Final positive-energy momentum distributions W v > = \<&^\^S?(T) >| 2 after the interaction with two 
counterpropagating laser pulses (£ = 1, T = 150tt/uj). Due to its magnetic component the field transfers longitudinal momentum 
to the wave-packet. Shown are the results for four different frequencies corresponding to four peaks in the resonance spectrum 
in Fig. [3] Due to the symmetry of the spectra under momentum inversion p' z — > — p z , we only show the positive half of them. 



oscillation pattern is strongly modified by inclusion of 
the laser magnetic field, which is seen most clearly in 
Fig. [2](b). The probability oscillates around a plateau 
value with a frequency five times smaller than the Rabi 
frequency Q in OEF. Note that the pair production rate 
at small times t <C scales as Q 2 [5[ and is reduced due 
to the magnetic field effect by 1-2 orders of magnitude. 
Taking for each frequency the maximum of the produc- 
tion probability as a function of laser pulse duration, one 
arrives at the spectrum in Fig. 03 In that way the nor- 
mally superimposed pulse length-dependent oscillations 
of the probability are omitted, and the resonances are 
clearly visible. In the OEF case, only odd-n resonances 
occur for p = 0, as the even ones are forbidden by a 
charge-conjugation related selection rule [l3j]. Due to the 
non-zero wave-packet width Ap, even-n resonances do 
occur in Fig. O but they are suppressed as compared 
to the odd ones. The picture changes significantly when 
real laser fields are applied: The height of the probability 
spectrum is reduced by approximately one order of mag- 
nitude, the resonances are shifted, several new resonances 
occur, and the resonance lines are split. 
In order to explain these modifications, we examined 
the corresponding momentum distributions of the cre- 
ated electron as shown in Fig. [4] for the peaks labeled 
by (3-1), (4-1), (2-2) and (3-2). In contrast to the OEF 
case, the photons in the CLP carry momentum along 
the beam axis, which is transferred to the electron wave- 
packet upon absorption. Only the transversal momentum 
components are conserved here. By energy- momentum 
conservation, a number of n + (n_) photons absorbed 
from the beam propagating to the right (left) determine 
the laser-dressed final particle 4-quasi-momenta: 

q' = n + /c + + ri-k- — q, (1) 

where q is the electron initial 4-quasi-momentum, k± = 
(cj,0,0,±cj) the photon 4- momenta for the right/left 



propagating waves. Our numerical calculations of the 
final momentum distribution after the laser fields have 
passed in Fig. [H confirm Eq. (pQ): the mainly contributing 
region for each peak corresponds to p' z = (n+ — n_)o;, 
for example p' z ~ 2uo at uo = 0.735m. The latter means 
that the final dressed momenta q' z do not differ essentially 
from the momenta outside the laser field p' z . In order to 
determine the resonance frequencies, we assume that the 
effective mass m* depends only on £ as in the OEF case 
(q 2 = ml) and that the initial quasi- momentum van- 
ishes (q = 0) because of the initial vanishing momentum 
p = 0. Then, the resonance frequencies read 

m* n + + n_ 
uo — — . (2) 

The main contribution to the momentum spectrum for 
the peak at uo = 1.1m in Fig.[3]comes from p' z = 0, when 
the number of absorbed photons from the left and right 
laser beam are the same. This peak belongs to the low- 
est possible photon number n+ = n_ = 1, resulting in 
m* = 1.11m. Then all peaks can be explained by Eq. ([2|) 
and correspond to the given labeling n_). The in- 
terchange of n + <-> n_ gives the same resonant frequency; 
we chose n + > n_ for the labeling. All resonances with 
n_ = 1 lie above uo = m*/2, according to Eq. ([2|). We 
found them up to a photon number of n + + 1 = 7. Their 
height decreases with increasing photon number. 
So far we have explained the overall structure of the res- 
onance spectrum in Fig. 0J There is however a substruc- 
ture inherent to all resonance peaks with n + ^ n_. The 
resonances are split into double peaks which do not occur 
for an OEF. The reason for this is again the spatial de- 
pendence of the laser fields so that the photons carry mo- 
mentum in propagation direction. The negative-energy 
electron can absorb n + photons from the left beam and 
n_ from the right, or vice versa, reaching two different 
final positive energy states with equal energy but oppo- 
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teraction [161 ]. 

In conclusion, we employed a numerical approach to in- 
vestigate e + e~ pair creation in counterpropagating laser 
fields. In comparison to the OEF case, we find character- 
istic modifications of the created particle spectra and the 
Rabi oscillation dynamics. The narrow peak splitting of 
the resonant pair production probability could serve as a 
sensitive probe of the quasi-energy band structure and, 
generally, of QED in superstrong spatially and tempo- 
rally inhomogeneous fields. In particular, the potential 
relevance of non-standard model field theories [17| could 
be tested. 

We acknowledge enlightening discussions with D. Bauer. 



FIG. 5: Coupling between the two upper states of the V- 
type system via Compton-scattering for the (2 — 1) peak, 
a) Population of the initially empty mirror state with p' z — 
— ^initial = — (solid line) co = 0.766m and £ = 0.5, result- 
ing in the Rabi frequency of = 0.018m due to Compton 
scattering; (dashed line) uo — 0.812m and £ = 1.0 resulting in 
Qc = 0.066m. In both cases Q <C He, where Q is the associ- 
ated Rabi frequency due to pair creation, b) £ dependence of 
the peak splitting, Aco ~ £ 2 . c) Level scheme. 

site momentum. Therefore, the pure two-level system 
of the OEF is transformed into a V-type three-level sys- 
tem (if n + 7^ ri-) in CLP, see Fig. 0(c). The positive- 
energy states are coupled via Compton scattering as in 
the quantum regime of free electron lasers. This can lead 
to splitting of these levels and, consequently, to splitting 
of the resonant transition from the^initial state, analo- 
gous to the Autler-Townes effect [14[. We investigated 
the splitting for the (2- 1) resonance. Fig. [5] (a) shows 
the Compton oscillation between the two positive-energy 
states for two different £ values. The Rabi frequency Qc 
due to Compton scattering increases by a factor of 3.6 
when £ is varied from 0.5 to 1. The splitting for var- 
ious £ is shown in Fig. 0(b) . A quadratic dependence 
Auo ~ £ 2 is found, leading to an increase of the splitting 
from £ = 0.5 to £ = 1.0 by a factor of 3.5. This indicates 
that the observed peak fine structure of the spectra is 
indeed an Autler-Townes-like effect. An equivalent ex- 
planation of the peak splitting can also be offered: The 
spatial periodicity of the field induces a band structure of 
the electron energies [HI which inhibits electron creation 
in certain energy regions. The splitting becomes larger 
with increasing £, following the energy gap behavior. 
Laser sources which could in principle be employed to 
enter the parameter range of interest are the XFEL fa- 
cilities presently under construction at SLAC and DESY 
(Hamburg, Germany). They are envisaged to provide co- 
herent synchrotron radiation with photon energies up to 
12 keV and could reach the £ ~ 1 regime if the laser beam 
is focused to sub-/im waist size [7]. Powerful hard x-ray 
sources are also expected via relativistic laser-plasma in- 
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